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CALCULATION  OF  IMPEDANCE  OF  A SHARP  PLASMA 
BOUNDARY  WITH  A MIXED-TYPE  ELECTRON  SCATTERING 
IN  ANOMALOUS  SKIN-EFFECT  CONDITIONS 

V. I. Miroshnichenko1,  V.M.Ostroushko2 

'institute  of  Applied  Physics  NASU,  ul.  Petropavlivska  58,  Sumy  40030,  Ukraine 
2National  Scientific  Center  "‘Kharkov  Institute  of  Physics  and  Technology” 
ul.  Academichna  1,  Kharkov  61108,  Ukraine.  E-mail:  ostroushko-v@kipt.kharkov.ua 

Some  aspects  of  calculation  of  impedance  change  caused  by  oblique  incidence  of 
electromagnetic  wave  on  a sharp  plasma  boundary  with  mixed-type  electron  scattering 
in  anomalous  skin  effect  (ASE)  conditions  are  considered.  Though  the  corresponding 
Riemann-Hilbert  boundary  problem  can  be  reduced  [1]  to  the  integral  equations  of  the 
Fredholm  type,  formal  limit  transition  to  extreme  ASE  causes  kernels  singularity.  The 
transformation  of  the  equations  to  ones  with  bounded  kernels  is  made  by  means  of 
partial  inversion  of  the  integral  operator. 

An  oblique  electromagnetic  wave  reflection  from  a sharp  plasma  boundary  with  a 
mixed-type  electron  scattering  can  be  described  with  the  functions  TT(kz)  and  Tt(k7), 
the  linear  combinations  of  the  Fourier  transforms  of  electric  field  components  in  plasma 
domain  (z>0).  They  should  be  analytical  in  the  half-plane  lm(kz)<0,  together  with  the 
functions  <J>>.(-kz)  and  <Eh(-kz),  and  should  meet  two  functional  equations  [2], 

[l-Q2q,(Pkxz)][^(kz)+p^,(-kz)]=  <Dx(kz), 

[l-kxz2-n2qT(pkxz)]['FT(kz)-pTt(-kz)]=(DT(kz),  (1) 

and  the  equalities  Tfr(±ikx)=±ivFT(±ikx).  Here  p=vFro[c  (m+iv)]‘],  Q=coe[to(©+vv)]_1/2, 
kxz=(kx2+kz2)1/2,  c is  the  speed  of  light,  to  is  the  frequency,  ooe  is  the  plasma  frequency,  v 
is  the  collision  frequency,  vF  is  the  velocity  of  electrons  at  the  Fermi  surface  (spherical), 
kx  is  dimensionless  (in  units  oo/c)  transverse  wave  number  (sine  of  wave  incidence 
angle),  pe(0,l)  is  the  fraction  of  electrons  with  specular  boundary  scattering, 

q^(x)=3x_2{(2xr1ln[(l+x)/(l-x)]-l},qx(x)=3(2x2)-,{l-(2Xr1(l-x>[(l+X)/(l-X)]}- 

The  problem  is  considered  under  the  approximation  P«l,  A»l,  where  A=(PQ)2/3,  so, 
the  distance  traveled  by  an  electron  during  the  oscillation  period  is  large  compared  with 
the  effective  field  penetration  depth. 

Introducing  new  variables  and  functions,  one  can  write  (1)  in  the  form 

X^.(u)+pQxx(u)Xx.(-u)= Y x.(-u),  XT(u)  -pQTx(u)XT(-u)=YT(-u),  (2) 

where  Xx(u)=  TfrOfr1  ^x((ipu)_l)Qx+(u),  Xx  (u)='Tt(-ctu)Qt+(u),  Cr=-ic3Ap_1, 
c3=exp(i7r/3).  Here  and  further,  a function  with  the  upper  index  + should  be  analytic  in 
the  half-plane  lm(u)>0,  and  it  should  be  built  for  the  corresponding  function  without 
upper  index  in  accordance  with  the  identity  Q+(u)Q+(— u)=Q(u).  A corresponding 
function  with  the  upper  index  x is  determined  by  the  identity  Qx(u)=Q+(u)[Q+(-u)]  ^. 
The  functions  Qx(u)  and  QT(u)  correspond  to  the  first  factors  in  the  (1)  left  hand  sides 
multiplied  with  the  constants.  The  identity  for  Q+(u)  remains  its  sign  indefinite,  and 
therefore  the  signs  are  chosen  to  provide  Q^+(ioo)=l  and  limu^ia)[QI+(u)/u]=cT. 
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By  using  (here  and  further)  the  tilde  symbol  for  the  Fourier  transforms  over  all  real  u 
axes,  from  the  first  of  (2)  one  can  obtain  the  following  equation: 

oo 

xAc)*p\diQl(C  + ^)xM)+pQ’A()=^  (t>»).  (3) 

o 

(3)  can  be  solved  by  iterations  [3],  and  some  values,  in  particular,  the  quantity  F*.= 
hmu_io0{iu[Q>L+(u)  -XAu)]},  can  be  efficiently  calculated.  The  function  XT(u)  can  be 
written  as  linear  combination  of  the  functions  X(±p;u)  determined  by  the  equations 

X(p;;u)  -p/Qtx(u)X(p/;  -u)  =Y(p/;  -u)  (4) 

(for  p =p  and  for  p/=-p),  and  with  the  requirements  that  the  functions  X(p  ;u)  and 
Y(p,u)  are  analytic  and  bounded  in  the  half-plane  lm(u)>0,  and  the  condition  X(p;; 
(oo)=l  is  satisfied.  From  (4),  an  equation  similar  to  (3)  can  be  obtained  and  the  quantity 
4^1=1  liml,_>iCO[uX(-p;u)  - cT-1  Qx+(u)j  can  be  calculated,  however  to  calculate  some 
other  quantities  one  has  to  transform  the  equation.  Taking  a> 0,  for  the  Fourier 
transforms  of  the  functions  x(a,p';u)  = X(p/;u)(ia+u)_',  g(a  ;u)=Qx>'(u)(ifl-u)(ia-i-u)“l, 
h(  a ,p  ;u)  = Y(p7;i  a )(in  +u)_1,  one  can  obtain  an  integral  equation: 

co 

x{a,p'-,£)=  p'\d%g{a\<Z + Z)x{a,p'\$)+h(a,p'\£)  (£>0).  (5) 

o 

The  function  for  1 «<^« A is  close  to  {nQ~\  so,  in  the  limit  of  A-*co,  (5)  is 

not  a Fredholm-type  equation.  For  p;=p  and  bounded  A,  one  can  take  a= 1 and  write  the 

A 

equation  x(l,  p;£)  = p + %)]  1 x(l,  /?;£)  + //(£”),  where  the  term  //(^)  involves 

i 

the  integrals  with  the  unknown  function  x(l ,p;£)  and  kernels  bounded  with  some 
functions  ([(C+^r2ln(^)+A-2(C+^)]  for  {t;e(l,A)£e(l,A)},  A(C+^)~2  for 

max(^,£,)>A,  and  max(l,^,£,)  1 for  other  domains.  In  brackets,  here  and  further  the 
factors  independent  on  and  A are  omitted.  Considering  the  function  H( £)  as 
known,  one  can  build  the  solution  of  the  last  equation  applying  the  second  of  two 
methods  used  in  [4].  The  kernel  of  the  corresponding  regularization  operator  is  bounded 
with  the  product  min(£',£)~*'’  max(^,£)'''’~1  (where  Kp=7if 'arcsin(p)),  and  the  kernel  of 
the  obtained  equation,  for  £e(l,A),  is  bounded  with  the  functions 
C, " " + A*-1’  s’ , and  d'+hr £ ^ intervals  ^e(O.l),  ^e(l,A),  and  ^>A, 

respectively.  The  introduction  of  the  variable  C,  and  the  function  x(^),  so  that  the 
products  of  derivative  dC,  jd£  with  some  factors  (1  for  i^g(0,1),  g2~Kl’  for  <^e(l.  A"'' ), 
A{C,  / A)'  k,‘  for  C,e(A*r  ,A),  C,1  / A for  C,  > A)  leads  to  the  bounded  limits  for  fixed  g 
if  A— mx>.  The  value  x(<^)/x(l, is  equal  to  1,  Kp , and  ^A~Kp  in  the  intervals 

^e(0,l),  £e(l,A),  and  Q> A,  respectively.  All  this  yields  an  equation  with  the  kernel 
bounded  even  for  A— >oo. 

In  the  case  p=-p,  to  keep  the  sign  in  the  estimation  of  the  product  p g(a;^)  for 
1«<£«A,  it  is  reasonable  to  take  positive  5j~A_1  and  consider  (5)  for  a = SX . Denoting 
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ca=C32A 5o=-(cCT1/2,  Ci=-tX(p;i5o)<50'f'’ , Co=-lX(-p;f8o) SQ  Kp , and  using  the  Euler 
gamma-function,  for  1 «Ct«A,  one  can  obtain  x(l,p;C)r(icp)»  , 

x(s„-p-x) r(i  — k p ) ~ C0^  Kp . In  the  limit  of  A— mx>,  the  quantities  Ci  and  Co  tend  to 

finite  real  values,  which  can  be  effectively  calculated  from  the  solutions  of  the  equations 
with  bounded  kernels  obtained  from  (5). 

Suppose  that  the  functions  XCJ(p/;u)  and  YCT(p/;u)  are  analytic  and  bounded  in  the  half- 
plane lm(u)>0  and  obey  the  functional  equation 

XCT(p/;u)  -p  QxX(ca/u)XCT(p/;  -u)  =Y0(p/;  -u)  (6) 

and  the  condition  XCT(p/;  (oo)=l.  Then  the  identity  X(p/;u)=  X(p/;0)Xa(p/;cCT/u)  takes 
place.  Applying  the  limit  transition  to  the  explicit  solution  of  the  equation,  in  which  the 
function  QTx(u)  is  approximated  with  the  relation  of  the  polynomials,  one  can  obtain  the 
equality  X(p;0)X(-p;0)=l.  From  (6),  an  equation  similar  to  (3)  can  be  obtained  and  the 
quantity  FCT=i  limu.„>iaj{u2(<3/ou)  ln[X  a(-p;u)/QU(cc/u)]}  can  be  calculated.  Besides,  (6) 
can  be  treated  similarly  to  (4),  and  can  equation  similar  to  (5),  but  with  another  function 

instead  of  g{a\£  + £),  may  be  obtained.  For  the  quantities  C0i=-iX0(-p;18o) Sqp  , 
C0o=-iX0(p;i8o)  S0  Kp  (which  one  can  effectively  calculate  considering  now  the  cases 
p=-p,  a- 1 and  p=p,  a-S]  and  evaluating  the  solutions  of  the  corresponding  integral 

equations  for  large  Q,  and  for  the  quantity  Fx  =X(-p,0)(c3A)  '' , the  identities 
CoCffr'=Ft  =CCToCr'  take  place. 

Under  the  assumption  (A»l,  p«l},  for  dimensionless  impedance  Z one  can  obtain 
approximate  identity  P(AZr'*vFTiexp(i7r/3)+kx2p2A"1"2'r'’  exp(-2i7iKp/3)  Fz,  in  which  the 
quantities  'Fti  and  Fz=(Fv-FCT)Fx2  in  the  limit  (P~ >0,  A— >co}  turn  positive.  The  complex 
factors  in  the  expression  for  the  impedance  are  obtained  explicitly  due  to  introduction  of 
convenient  variables,  for  which  the  corresponding  functions  in  the  functional  equations 
in  the  extreme  ASE  limit  become  positive  for  real  values  of  arguments. 

The  quantities  Ci,  Co,  and  Tfi  can  also  be  calculated  with  the  method  of  [5]. 
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